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1.0 


INTRODUCTION AND SUMMARY 


This report documents the results of several Shuttle/GPS studies 
conducted by Axiomatix under Contract NAS 9-15387B for NASA's Johnson 
Space Center. These studies augment the studies conducted by Axiomatix 
for NASA under the companion contract NAS 9-15387C. Several of the 
investigations reported here, such as the GPS system self-test, the 
ground transmitter, and the ionospheric modeling investigations, grew 
out of needs which were developed In the Shuttle/GPS Panel meetings that 
Axiomatix participated in. The effects of oscillator noise and the tran- 
sient maneuver analyses were conducted as part of the effort to determine 
the GPS receiver performance in the Shuttle environment. These studies 
and investigations are summarized briefly below. 

The investigation of the effects of oscillator noise on the GPS 
receiver developed comprehensive models of the oscillator phase noise 
process and the phase lock loop response to phase noise. It is shown 
that, for a reasonably low phase noise oscillator, the degradation to 
Shuttle GPS receiver performance is negligible. 

A Shuttle/GPS receiver self-test mode is investigated. This 
investigation was motivated by the desire to utilize the one good GPS 
receiver in the doubly redundant system in the event of one failure. 

A viable technique is developed that is capable of testing the GPS/ 
Shuttle data interfaces, as well as the receiver itself. 

Since the Shuttle will utilize GPS before the full constellation 
of 24 satellites is deployed, it may be desirable to use a GPS ground 
transmitter for Shuttle landing operations. Such a ground transmitter 
is described and a cost estimate made. 

The GPS system utilizes two frequencies, 11 and L2, to correct 
for ionospheric delay. An investigation is performed to determine if 
the ionospheric delay could be modeled so that dependency on L2 might 
be eliminated. The conclusion is reached that this is not practical and 
is a poor receiver complexity trade-off. 

Since the Shuttle fires thrusters to initiate maneuvers while in 
orbit, transients of range acceleration or jerk occur. Such transients 
can cause narrowband receivers to lose lock. Thus, the performance of 
a Shuttle/GPS receiver with typical Shuttle transient maneuvers has been 
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analyzed. Worst-case receiver bandwidths were used and it is shown that 
the receiver tracking performance is negligibly affected. 
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2.0 EFFECTS OF OSCILLATOR NOISE ON SHUTTLE GPS 

RECEIVER PERFORMANCE 

Introduction 

Oscillator phase noise, generally referred to as 1/f noise, 
degrades the doppler measurement capability of the GPS receiver and, 
if severe enough, can raise the receiver threshold. This latter effect 
can cause the receiver to lose lock prematurely, lose data, or take 
longer for acquisition of the PN signal. The extent of receiver degra- 
dation is dependent on the receiver carrier tracking loop bandwidth and 
the phase noise spectral density of the carrier. Thus, it is important 
to analyze the carrier (Costas) tracking loop phase error caused by 1/f 
noise for nominal Shuttle GPS receiver bandwidths and a typical "good" 
oscillator phase noise spectral density. The following analysis has thus 
been performed, and the results show that 1/f noise is not. expected to 
degrade receiver tracking or data performance. 

Analysis 

The model of the phase locked carrier tracking loop used for the 
analysis is given in Figure 1, where 

S(t) = v^A sin [co 0 t+e 0 (t)] 
and A = amplitude of input signal 

* o n (t) = 1n P ut Phase be constant or may represent doppler 

1 or frequency offset) 

r(t) = reference voltage = IT cos [iDgt+ 6g( t)] 

K, = voltage gain 

> Sg(t) = P^ ase estimat e 

e ( t ) = dynamic phase = K^y(t) r(t) (at baseband) 

= « d A sin <t( t) 

K d = detector gain 

» 4(t) = o Q ( t) - 0 Q ( t) = phase error 

The VCO operation is given by 

se Q = Ky[z( t) a. v. ] (K v in rad/sec/volt) . 

9 

Assume noise = 0 and acquisition voltage (a.v.) = 0, so that 


9 
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e n~ $ 


K v K f(s) 

_v z (t) = — e(s) 


= A K 1 K d K v 


sin 4 >(s) . 


Let K 4 K, K d K y 


loop gain and thus $(t) = 9 q( t) - A K "f ( s ) sin $(s) 


iM * d 9 o (t) AK '* 


dt 


dt 


f ( t - t) sin <t> ( t ) dt . 


( 1 ) 


Equation (1) is the nonlinear integral equation that represents operation 
of the phase lock loop. Its equivalent mathematical model is shown in 
Figure 2. 



Figure 2. Mathematical Model of Nonlinear Phase Lock Loop 


Assume the phase error <p is very small; then sin <j> = <p and the 
loop will be linearized and will take the form shown in Figure 3. 



Figure 3. Linearized Tracking Loop 




The linearized closed-loop transfer function H (s) * ®q( s )/ e Q ( s) from 
Figure 3 Is 

( 0 . § ) hJUtel * § 
ie 0 V s V 

L e t Q » - - i then Ge Q - Ge Q 3 e Q , 6^(1 + G) = G6 q, and 


H(s) 


8 fO a _G 


A K f ( si 


1 + G s + A K f(s) ' 


( 2 ) 


The relationship between 4>(s) and e Q (s) is 


e 0 ” ^G = $ 


*(1 + G) - e Q 


4 = 


Jj 3 

+ GJ 3 0 ’ 


and from (2), 


*(s) 3 [1 - H(s)3 e Q (s) 


(3) 


The equivalent filter is shown in Figure 4. 



Figure 4. Equivalent Noise Filter 
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Assume the power spectral density of the error process <j>(t) to 
be S^(f), and S 0() (f) Is the power spectral density of the error process 
e Q (t). Then, 

S.(f) * S (f) |l-H(ju.)l 2 . (4) 

♦ °0 

5 If we assume the process ®qU) has a null mean value, then 

0 o (t) - e Q (t) 8 $(t) 8 0 

> and the variance of the error phase process 4 >(t) Is 

o 2 8 f S 6 (f) 1 1 - H(jai) | 2 df . (5) 

♦ e 0 

It Is required to calculate o 2 using (5) under the following 
assumptions: The loop filter is as shown in Figure 5. 


K 
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Figure 5. Loop Filter 
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The loop filter transfer function F(jw) can be calculated as follow? [1,2]: 

U ? (ju>) l + ju)T ? 

F(jw) 3 7 - ^r a -G tt 7 — — r . 

U-j IJwJ 1 + ju) T-j 

where t 2 = R 2 C 

3 + GR 1 + R 2 ) C . 

Assume G is relatively large; then 


GR 1 C . 


l+jwTj, 1 + ju T 2 1 + jtO T 2 

= 6 1 + jw G R|C : ' ju» R 1 C = jw rj“ 


( 6 ) 


where t 2 = R 2 C and t 1 = R^C 


by [1]: 


The natural angular frequency of the loop u> n is always given 


AK 

T 1 


(7) 


and the damping constant c by 


c = \ t 2 


( 8 ) 


This implies 


2f, = A K 1 2 / t ^ 


(9) 


Substituting F(j w ) of (6) into (2) gives 

A K ( 1 + j w 1 2 ) 

H(J “’ = j4>tjy+A'kir+j7r 2 7 ■ 


AK/ Tl + jo) AK t 2 / t-j 
♦ AK/ T ^ + juiA K T 2 / t ^ 


Using (7) and (9), we get: 


9 
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H(j«) 


+ 2 ju> X, UJp 

~1 — i — ; 

-u + U) n + 2 j 10 £ U)^ 


and 


1 - H( joi) = 


2 . 9 , r 2 2 , , r 

+ C J U) £ " W “ U)^ - C J 0) ? 0)^ 

T 


0 ^ + 2 jiu c 0) n - 0) 


2 9 ^ r 2 

0 - i Jw c 0> n - U) n 


|1 - H(jo>) r = (l “ H(ju)) (1 - H(ju)) 


0) U) 

^ x ^ g 

0) " U)^ - 2 JCO £ CJ^ CO - U) + 2 j (a) £ < 1 ^ 


^ 

(co 2 - co n 2 ) 2 + (2w c 0J n ) 2 


4. 4 9 2 2. , 2 2 2 * 

CO + 0)^ “ COJ + HCO 0) £ 


n 


n 


Assuming c = . we get 

Ji 


n-H(j U )r 


4 . 4 

0 + U) 


( 10 ) 


It is reasonable to assume that the oscillator phase noise spec- 
tral density S 0 g(f) has the shape given in Figure 6 and denoted by L(f). 
From Figure 6 , we get 


s e < f > ' ^ 


11.7 


0 < f .c 4.3 Hz 


_ 10 ‘ 


•13.0 


> 4.3 i f < 30 Hz 


10 


-14.45 


(ID 


f > 30 Hz 



Jm 
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Substituting values of (10) and (11) In (5) gives 


s 2 


4-3.Q-11.7 


u) 4 

+ CD 4 

n 


df + 


30 10 -13.0 
4.3 f 


w 


df 


w + U) 

n 


f B IF , 0 -l«.«5 J_ ' 
J 30 w + u 


Ui 1 + i 2 + l3 ]. 


Equation (5), the tracking error variance, will be evaluated for a narrow- 
band tracking loop, representative of a code tracking loop, where 
u n =2 rad/sec, and for a wider tracking loop, w n *30 rad/sec, representa- 
tive of the carrier tracking loop. In both cases, the receiver IF 
bandwidth is taken as 30 Hz. 

The values of integrals Ij, Ig, and I 3 are calculated using tables 
of integrals [3] as follows: 


= 10‘ 11,7 ( 2 n ) 2 


r4.3 


a) 


o- ^ _ 7.87698 

-17T 2 ” df - — rr- 

<i>+ w 10 


f8.6iT 


U) 


16 + w 


du> , 


and, from [3, p. 62], 


a + bx 


dx = — - — tan’^ x 2 /b/a , 


2/iT 


so that 


t - 7.87698 f 1 f -1 ,2 /a : 

*1 ' 1Q 11 [2x4 tan ^ /4 l 


8.6tt 


= 1.54x10' 


■11 


Similarly, 


\ 2 - 10 


■13.0 


2ttx30 3 


u> 


2ttx 4. 3 16 +w 


du> , 


and from [3], 


dx 


a + bx 


4b 


In (a + bx 4 ) , 
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l 


* 


# 


>2 ‘ 


i a^I 3 • D 1 i f i c a. ^ \ 

10 ^ in (16 + to ) 


-i 60 tt 


1 .94 x TO” 13 . 


8.6tt 


Likewise, 


10-' 4 « £ f" V-T 2ir df 

" J30 w +uT 


0.05647 x 10 


■14 


f400ir 4 
■'60 tt 16 + oj 


^ dto . 


and, from [3], 


r f — ~ — r dx = lfl±^ 
b * a + bx 4 b ' a + 


bx 4 -_a dx , I 

4 ax b 
a + bx^ D 


dx- f — 
J a + 


dx 


a + bx 


where 


1 


' a + bx 


dx = 


Va7b 

4a/2~ 


y t -M ♦ 2 

x 2 -Vi7bx/2 + /57F Wb- x 2 


)] 


Thus, 


I 3 = 0.05647 x 10 


-14 


0) - 


2/2 


fin ♦ 2 tan' 

to - 2/2 0) + 4 


1 / 2/2 U)\ 

14-4 


n400ir 
60 it 


= 6x 10 


-14 


Upon combining the three integrals, the rms tracking error is 
found to be 

o = 5.7x lo~ 6 rad . 

4> 

Proceeding in a similar manner, the rms error for the to^ = 30 rad/sec case 
is found to be 1.3xlo" b rad. 

It is important to note that these errors are calculated for a 
5 MHz oscillator. 
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Since the GPS LI frequency is 1.575 GHz, the oscillator is 

followed by a multiplication function so that the equivalent receiver 

g 

oscillator has the frequency 1.575 x 10 Hz, as shown in Figure 7, 


y 


h(f) 


vco 


Operating 
frequency 
is 5 MHz 


Figure 7. Equivalent VCO Operating at 
1. 575x1 0 9 Hz 


where h(f) = — — - x | — ■ = 315 
5x10° 

S y (f) = |h(f)| 2 S x (f) = (315) 2 S x (f), 

so a correction on error variance has to be done for both cases as 
follows: 

Case 1 u> n = 2 rad /sec, B Ifr = 200 Hz/sec 

o 2 = 32.420793 xl O' 1 2 x (31 5) 2 

a = 0.0017935 rad = 0.10256 deg. 

<P 

Case 2 oj_ = 30 rad/sec, B ir = 200 Hz/sec 

n Ir 


o 2 = 160.82877 x 10 -14 x (315) 2 

= 0.00039947 rad = 0.02289 deg . 


Thus, it is seen that the phase error due to oscillator noise is negligibl 
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3,0 GPS NAVIGATION SYSTEM SELF-TEST 

System self-test is a concept that enables the GPS Navigation 
System to validate its own performance in real time during a mission. 

This concept should not be confused with certain avionics Built-In Test 
E uipment (BITE) which is usually used to isolate processor faults to the 
irtdule or board level. Also, BITE is sometimes a software function used 
•:o verify the correct operation of receiver software and firmware. An 
example of such a BITE function for a GPS receiver is given in Table 1. 


Table 1. Sequential Receiver Built-In Test Concepts 
for Software/Firmware 


e Power Fail/Restart Interrupt 

- Performs critical data memory validity check on restart 

- Calls system initialization on restart 

- Flag made available to distinguish fail from restart 

• Computation Error Interrupts 

- Floating point microcode 

- FORTRAN library errors 

- Errors tabulated and computation continued or process restarted, 
depending on error 

• Time-Out Error 

- Used fc. stack overflow, jump to restart 

• Self-Test Background Task 

- AFI monitoring 

- ROM checksum 

- RAM read/store/restore/compare tests 

- Critical data memory checksum 

• Each I/O Ta'.k Provides Bit Functions 

- Telemetry 

- Receiver 


Since the baseline Shuttle GPS Navigation System is a two-string 
system, when thi two GPS navigation solutions do not agree due to a mal- 
function ir one of the strings, both are discarded and the Shuttle 
utili.fs the baseline navigation system (i.e., TDRSS doppler, or TACAN). 
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An effective self-test feature would disclose which GPS receiver was 
still providing good navigation data and enable this data to be utilized 
by the Shuttle, if so desired. 

The desirable features of a self-test function include the abil- 
ity to provide a system end-to-end test which provides an absolute indi- 
cation of the system health of the entire GPS navigation chain. This 
indication should be a continuous, real time, in-flight health evaluation. 
However, it should also be capable of being used as a maintenance tool 
to check out the receiver performance prior to installation in the Shut- 
tle. Each GPS receiver should have its own self-tester, independently 
operated. The self-test should provide a verification of receiver "front 
end" health, including tracking and data functions, and pseudo range 
and doppler measurement functions. It is also possible to consider 
verification of receiver processor and memory health. This can be 
accomplished by a "test state vector" and "test bite" indicator which 
are compared in the Shuttle GPC with the correct design values. 

A block diagram of a possible self-test implementation for a 
sequential GPS receiver is shown in Figure 8. The essence of this sys- 
tem is the internal generation of a GPS-like signal that the receiver 
processes. In processing this signal, the receiver must exercise, and 
thereby check out, all the functions required for normal GPS satellite 
signal reception. The basics of the test GPS signal generation are 
shown in the block diagram in Figure 9. As can be seen, the self-tester 
employs a GPS signal generator that modulates an LI carrier with a P 
code and appropriate data stream. This signal is injected, at an appro- 
priate RF level, as close as possible to the system antenna. The two 
most logical places to inject the test signal are at the preamp input 
and at the receiver input. The advantage of injecting the signal at 
the preamp input is that the entire RF circuit is tested; however, this 
is accomplished at the expense of an additional coax cable which must 
be run from the receiver to the preamp. Even though this cable can be 
a small -diameter, relatively lossy cable, the running of an additional 
cable to the preamp can present problems due to severe space limitations. 
Also, this signal injection at the preamp input will cause a slight 
increase in the preamp noise temperature. 






RECEIVEO 
GPS SIGNAL 



Figur? 9. Self-Tester Functional Diagram 
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Injection of the test signal at the receiver input avoids the 
problem of running an extra coax cable to the vicinity of the antenna. 

It also avoids the additional increase noise temperature; however, It 
does not provide for checking the preamp. Thus, if a preamp falls, the 
self-test may still indicate satisfactory operation of the system. This 
is not as serious a problem as It might seem since preamplifier failures 
are generally "hard failures" which are detected by other receiver health 
indicators, such as receiver AGC. Thus, it would seem that the receiver 
front end is a reasonable point to inject the self-test signal. 

The self-test function must be automatically time-multiplexed 
with the normal GPS signal measurement function. This concept is illus- 
trated in Figure 10 for a sequential receiver. Normally, the sequential 
receiver sequences from satellite number 1 slot to satellite number 2 
slot and so on until four measurements have been completed. The sequence 
is then back to satellite number 1, with appropriate substitution of new 
satellites as they become visib’e. A typical frame length is no more 
than 24 seconds long. The test function is multiplexed by creating a 
fifth slot in the frame for testing Since it is not necessary to verify 
system health as frequently as every 30 seconds, the fifth slot can be 
multiplexed itself with other processes such as search and acquisition 
of additional satellites. Thus, a typical sequence for slot 5 might 
be two frames of search for new satellites, followed by a test frame. 

This would provide health verification every 90 seconds or less. Slot 5 
could easily be software-configurable. 

A simplified flow diagram of the self-test process is shown in 
Figure 11. The test signal that is generated is designed such that the 
transmitted code phase and carrier doppler are preprogrammed to corre- 
spond to a desired range and range rate. Furthermore, the data message 
on the test signal would contain the ephemerides for a favorable geomet- 
rical position of the "test satellite." If it is desired to just verify 
receiver RF and baseband processing, only the pseudo range and delta range 
words and the data word need to be compared with the stored "test" 
words. A more complete test, however, of the receiver is to have the 
navigation processor portion of the receiver process the "test" pseudo 
range, delta range and ephemerides to calculate a "test" state vector. 
Since the comparison with the true value is performed in the Shuttle GPC, 
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Figure 11. Shuttle/GPS RPA Self-Test Flow Diagram 
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! this test also evaluates the GPS RPA/Shuttle data Interface. This Is 

the test that Is shown In Figure 11. 
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4.0 GPS GROUND TRANSMITTER DESIGN TO AUGMENT SHUTTLE NAVIGATION 

When the Shuttle enters the atmosphere for landing, there Is a 
possibility of significant enhancement of the landing navigation accu- 
racy with GPS by utilizing earth-based GPS transmitters. The earth-based 
GPS transmitters, hereafter referred to as ground SAT, offer a more 
favorable geometry for calculation of the GPS navigation solution. This 
is especially true for Shuttle operations prior to the fully deployed 
GPS 24-sa tel 1 1 te constellation. The ground SAT would be located within 
approximately several hundred miles of the Shuttle landing site. This 
section reports the results of Axiomatlx's investigation of the design, 
performance and cost of the ground SAT. 

The function of the ground SAT is to emulate a normal GPS satel- 
lite transmission. Of course, the data message regarding the ephemerides 
of the ground SAT is drastically different than that of the normal GPS 
satellite; however, the navigation signal structure is exactly the same. 

A functional block diagram of a ground SAT is shown in Figure 12. The 
major difference between this diagram and an equivalent diagram for a 
normal GPS satellite is the GPS user receiver that is present in the 
ground SAT. This receiver is used to synchronize the ground SAT clock 
to the GPS system time. To do this, the user receiver receives the nor- 
mal GPS navigation signals and solves for time, one of the parameters 
available from GPS. The receiver need not be a very sophisticated GPS 
receiver in terms of capability. This is because it is stationary, i.e., 
no dynamics, and that it can use a long period of time to converge to a 
highly accurate solution. A good candidate receiver is a Manpack 
receiver of a Z set. 

The oscillator shown in Figure 12 is a rubidium oscillator. The 
requirement for a high-quality oscillator stems from the fact that it 
may not be practical to simultaneously transmit the emulated GPS navi- 
gation signal and receive the normal GPS satellite signals at the same 
location. A possible trade-off is to remotely locate the GPS user 

i 

receiver and antenna so that they would not be saturated by the emulated 
signal. A cable of known length would then be used to carry the time 
signal to the oscillator. This has the disadvantage of requiring two 
sites and the inherent reduction in site flexibility. Furthermore, 

with the approach wherein a stable oscillator is used to "coast" during j 

1 

j 

j 






Figure 12. GPS Ground Transmitter for Shuttle Landing - Ground SAT 










24 


the transmission interval, one antenna with a suitable RF switch can be 
used for both transmit and receive functions. This is indicated by the 
dotted lines in Figure 12. 

A summary of the specifications for a commercially available 

rubidium oscillator is given in Table 2. It can be seen that, with a 

-12 

short-term stability of 1x10 (t=100 sec) and a long-term stability of 

<lxl0’^/month, this oscillator will provide excellent performance for 
the ground SAT application. The price of the oscillator is approximately 
$7000. 

The ground SAT would be housed in a portable shelter such as that 
depicted in Figure 13. Figure 13 actually shows an existing ground trans 
mitter used I : the GPS Joint Program Office (JPO) for the Yuma GPS 
inverted range tests. A block diagram of this transmitter is shown in 
Figure 14. This block diagram contains all the elements shown in Fig- 
ure 12, with the exception of the rubidium oscillator. The cost of the 
JPO ground transmitter has been reported by JPO to be approximately 
$100,000, including nonrecurring costs. Axiomatix estimates that the 
cost of a ground SAT for use with the Shuttle landing would be approxi- 
mately $200,000. This would include the rubidum oscillator, the GPS 
user receiver, and redundant elements for critical signal paths. 

The link budge for a ground SAT is given in Table 3. This budget 
is c.'culated for a 1 watt transmitter and includes 0.3 dB of antenna 
feed line loss. A range of 700 miles was arbitrarily chosen. The actual 
range will be less. A link margin of approximately 7 dB is obtained, 
which is more than adequate. 


t 
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Table 2. Rubidium Oscillator Performance 



Model FRK-L 

Model FRK-H 

Output 

10 MM/ sine wave 0 5V mis into 50 ohms, floating ground 
OPIION 5 MHz. 1 MHz 

Signal to Noise (SSB 1 Hz BW) 

>120 dB at 100 Hz and >145 dB at 1000 Hz from earner 

Input Power 

1 3W at 24 VDC, 25°C ambient; 22 to 32 VDC ; peak current during warm-up. 1 8A 


Warm-up Characteristics 


Long Term Stability (Drift) 
Short Term Stability 


Trim Range 


I nvironmental Effects 


Voltage Variation 


• Operating Temperature 


Storage Temperature 
Magnetic Field 


Altitude 


Humidity 


Shock 


Vibration 


General Information 


Size 


Weight 


Electrical Protection 


< 10 minutes to reat h l \ 10 10 at 25°L ambient 


< 4 x 10"’ ’/month 

3 x 10"’ 1 t = 1 set 
1 xIO •’ t ^ 10 set 

1x10 14 t - lot) sec 


< 1 x 10*' '/month 

1 x 10"” t = i sec 
4x10 ” t 10 sec 
I x 10 ” t -- 100 sec 


2 x 10"« 


1x10 


< 1 


<6 x 10 '° from - 40°C to + 6S°C 


x 10"’ 710% 


< 1 x 10'”’ from -25°C to + 65°C 


OPTION <1 x 10 * from -54°Cto 
+ 65°C 


OPTION <4x10 10 from - 54°C 
+ 65°C 


S4°C to + 75°C 


< 4 x 10"' ’/Am*’ (3 x 10"”/0 1 millitesla) Optional shield available 
< I x 10"' Vmbar (sea level to 21,000m) 


9S% MIL T 5422E 

Mil STD «10( . Method Sib 2, Proi edure I 
MIL STD-810C, Method 514 2, Procedure I 


100mm x 99mm x 1 1 2mm (3 9 in x 3 9 in x 4 4 in ) 


1 3 kg (2 9 lbs ), with optional heat sink EE K-10. 1 55 kg (3 5 lbs ) 

An internal diode and fuse protects against reversed polarity connection 


Connectors 


Coaxial connector OSM 211. mates with OSM 501-1 or equivalent 
i ight push-on connector pins 

OPTION Winchester Connector PN SRE-20PJ, mates with SRE-205J 


Warranty 


1 year, lamp and resonance cell-5 years 


•Operating temperature i$ baseplate temperature with optional heatsink EEK-10 attached Highest operating 
temperature depends on heat transfer between heat sink and unit 


OMNIDIRECTIONAL ANTENNA 
AND GROUND PLANE LOCATED 
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Figure 13. GPS Inverted Range Ground Transmitter 




Figure 14. Ground Transmitter Block Diagram 
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Table 3. GPS Ground Transmitter/Shuttle Link Ll-P 


Transmitter Output Power at Antenna 
Antenna Gain 
Space Loss 

Polarization and Atmospheric Loss 

Shuttle Antenna Gain 

Circuit Loss 

Received Power 

System Noise Temperature 

Boltzman's Constant 

Noise Spectral Density 

C/N 0 

Required C/N Q (Theoretical) 
Implementation Loss 
Required C/N Q 
Link Margin 


.3 dBW (1 watt transmitter) 

- 3.0 dB 

-153.5 dB (R= 700 miles) 

- 0.5 dB 

- 4.0 dB 

- 1.0 dB 
-162.3 dBW 

27.5 dBK 
-228.6 dB-W/K/Hz 
-201.1 dB-W/Hz 

35.8 dB 

29.6 dB/Hz 
2.0 dB 

31.6 dB/Hz 
7.2 dB 
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5.0 EFFECTS OF IONOSPHERIC DELAY MODELING ON GPS RECEIVER DESIGN 

Radio waves which pass through the earth's ionosphere travel more 
slowly than in free space due to group retardation along the ray path. 

This retardation is proportional to the total number of free electrons 
encountered and, as a consequence, varies with the angle of the particu- 
lar ray path, the time of day, time of year, and geographic location. 

The GPS navigation signals are likewise retarded by their propagation 
through the Ionosphere. This retardation manifests itself as an error 
in the measured pseudo range, which is the basic observable for computa- 
tion of the Shuttle's state vector with GPS. The approach taken to solve 
this problem has been to broadcast a navigation signal on two frequencies, 
LI and L2, with the user receiver measuring pseudo range on both frequen- 
cies. Since the delay is proportional to the inverse of carrier frequency, 
the user receiver is able to calculate the absolute delay and thus accu- 
rately determine the pseudo range. This section examines the feasibility 
of accurately modeling the absolute ionospheric delay, tnus eliminating 
the need for two frequencies. The impact on receiver design is examined 
along with the software impact to determine if there is a possible cost 
savings for the Shuttle GPS navigation system. 

An RF signal propagating through the ionosphere is delayed by 

K 

t = ~2 

Y 

where K is a parameter proportional to t;tal electron content along the 
measurement path and F is the signal frequency. If we consider two fre- 
quencies, F l -| and F^, the delay on LI is 

. . K 



and the delay on L2 is 

K 

t L 2 = 7T 
h L2 

The two-frequency receiver actually measures the differential delay between 
LI and L2, at. This can be written as 
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Solving for K, 



This last equation demonstrates the principle of the two-frequency GPS 
receiver. However, the immediate interest is to determine whether the 
ionospheric delay can be modeled accurately enough to eliminate one of 
the frequencies. 

The first ionospheric delay model studied is called the "slab 
thickness" model. It relates the delay to fOF2, the critical frequency 
of the F2 region of the ionosphere. The parameter fOF2 has been regularly 
measured for over 40 years by ionospheric sounders, called ionosondes. 
fOF2 is proportional to the electron density at the maximum of the F2 
region, called N_ av . A second parameter which has been measured for many 
years is MUF-F2, the maximum usable frequency of the F2 region, which has 
been shown to be closely related to the height of the maximum density of 
the F2 region. The ionospheric Total Electron Content (TEC), propor- 
tional to time delay, has been measured at only a few scattered locations 
via Faraday rotation effects. Klobuchar and Hawkins, however, have used 
available TEC data t-xen at Hamilton, Massachusetts, to form the ratio 
of TEC divided by N . which is called equivalent slab thickness. N 
values were obtained from fOF2 measurements made at Wallops Island, Vir- 
ginia. This slab thickness ratio of two experimentally measured quanti- 
ties is a first-order scale height or thickness parameter of the F2 iono- 
spheric region. Thus, if a value of fOF2 is known, the ionospheric time 
delay can be obtained from 


AT 


1 .67 x 10 “ 6 t 


(seconds) 

r c 
■ 1 1 


where fOF2 is in MHz, x comes from the slab thickness model and is 
expressed in km, and F^j is the GPS LI carrier frequency. 

A model for t was derived from the existing measured data. 


? 



A model value of t is determined from 


t = t(h. M, F10.7) = [M(month)* F 10.7 + B(M)] *C(h, M) 

where h is the hour of the day, M is the month of the year, and 10.7 is 
the 10.7 cm solar flux coefficient for a linear fit. These model coef- 
ficients are tabulated in Table 4. 

The error in this slab thickness model is shown plotted in Figure 15. 
The error is greatest, approximately 10%, during the night hours since the 
total ionospheric delay is least at these times. Conversely, the error 
is smallest, 5%, during midday when, due to the additional ionization 
caused by the sun, the ionospheric delay is greatest. Thus, it might be 
concluded that this is an excellent model; however, the serious drawbacks 
for this model are as follows: the model is appropriate only for a 

restricted geographic locale and, since approximately 4 kilobits of com- 
puter storage are required for the model, excessive computer storage 
would be required for a worldwide model. Furthermore, the research to 
extend the model to a worldwide basis has yet to be performed. 

Several models of ionospheric delay have been developed by The 
Johns Hopkins University Applied Physics Laboratory. An example of these 
models is an algorithm in which the ionospheric group delay is obtained 
via ray-tracing based on a priori estimates of the electron density dis- 
tribution. The ionospheric group time delay is given to the first order 


as 


where 


t(g, t) = 

F L1 


N(r, t) ds 


G(t) 


t(G, t) = ionospheric group time delay of the ranging signal 
between satellite and navigator along the geometric 
path G(t) at time t, ns 

f = carrier frequency, Hz 

N(r, t) = electron density at r and t, electrons/m 
r * position vector to a point on G(t) 
ds = differential element of arc along G(t), m 


A vertical electron density profile is used to characterize the 
electron density distribution. The model for this is shown in Figure 16. 
This model is characterized by 11 parameters, each of which can be a 
function of geographic location, time, solar activity, etc. The 11 
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Figure 16. Eleven-Parameter Vertical Electron 
Density Profile 
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parameters are determined In the following manner: 

(1) Density at 5000 km— an analytical function of the month of 
tha year and the gyrofrequency. 

(2) Density at 1000 km— an analytical function of magnetic dip 
and solar zenith angle. 

F2 region 

(3) Maximum electron density— obtained as a function of geographic 
location, unlversa 1 time, month, and solar index from praolctlons of f0F2 
available from N0AA. 

(4) Height of the maximum— obtained as a function of geographic 
location, universal time, month, and solar index from predictions of 

M{ 3000) F2 available from N0AA. 

(5) Semi thickness— obtained as a function of geographic location, 
universal time, month, and solar index from predictions of hmF2/ymF2 
available from N0AA. 

FI region 

(6) Maximum electron density- -obtained as a function of solar 
index and solar zenith angle from an empirical formula. 

(7) Height of maximum— obtained as a function of solar zenith 
angle from an empirical formula. 

(8) Semi thickness— constant at 50 km. 

E region 

(9) Maximum electron density— obtained as a function of geo- 
graphic location, universal time, month, and solar index from predictions 
of fOE available from N0AA. 

(10) Height of maximum— constant at 130 km. 

(11) Semithickness— constant at 20 km. 

The residuals for this, along with two other APL algorithms, are 
shown plotted for three locations in Figure 17. It can be seen that the 
residuals are on the order of 10 nanoseconds, which is an excellent result. 
As far as applying this model to Space Shuttle GPS navigation, however, 
there is the problem of gathering all the pertinent parameters discussed 
above, for a wide geographic area and time, and storing the data in the 
GPS receiver computer. Thus, the two-frequency receiver hardware penalty 
must be examined in light of this software penalty. 

The major impact on GPS receiver design of two-frequency, L1/L2, 
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operation Is in the front-end down-converter and In the receiver 
frequency synthesizer. Substantial Increases In receiver complexity are 
avoided by time multiplexing the receiver RF, IF, and baseband processing 
between the LI and L2 frequencies. Figure 18 shows a typical RF front-end 
and down-converter design for a two-frequency receiver. The Increased 
complexity due to LI and L2 frequencies Is shown by the dotted lines. 

The frequency synthesizer for the same typical GPS receiver design is 
shown In Figure 19. Again, the increased complexity due tu L1/L2 opera- 
tion is indicated by the dotted lines. The added complexity tor L1/L2 
operation, relative to the overall receiver complexity, is slight. When 
compared with the research, software, and computer firmware required to 
implement a suitable ionospheric delay model, the conclusion is reached 
to retain the two-frequency approach for the Shuttle GPS receiver. 
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6.0 GPS RECEIVER TRACKING OF SHUTTLE TRANSIENT MANEUVERS 

The Shuttle GPS receiver tracking ops must be sufficiently 
narrowband so that the tracking error due to thermal noise is small. 
However, this situation generally presents the classic coherent receiver 
trade-off between loop noise performance and dynamic maneuver tracking 
capability. That is to say, a narrowband loop is desirable from a noise 
consideration, and a wideband loop is desirable from a dynamic tracking 
consideration. The analysis which follows examines the dynamic tracking 
performance for the GPS receiver carrier and code loops. Typical values 
of bandwidth are used and the dynamics used are typical Orbiter on-orbit 
maneuvers. 

The acceleration rate, or jerk, profiles for an Orbiter pitch 
maneuver start and an OMS burn are shown in Figure 20. These profiles 
have been modeled as pulses, as shown by the dotted line in Figure 20. 
Although this is a worst-case model, it is analytically tractable. 

The tracking loops for both the carrier loop and the code loop 
are taken as second-order phase lock loops. The loop transfer function 
for the second-order phase lock loop is given by 



where 

r = AK *2 Fq 
F q = t 2 /x 1 

C = loop damping factor = /r~ /2 
= loop natural frequency = /r r /t 2 
= two-sided noise bandwidth = ( r+ 1 ) / 2 t ^ 

The loop error transfer function is given by 



o hLz J h l = 

e(s) 



s» 


1- H(s) 




# 


Figure 20. Jerk Models for Receiver Transient Performance Evaluation 
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A pulse in jerk in the time domain is given by 
.2 (t-t ,) 3 

e(t) = e j- u(t) - e — g L u ( t-t-j ) 

where u(t) is the unit step function and 


where 


° - Vo 

a = acceleration rate or jerk in ft/ sec.' 
c * propagation velocity in ft/sec. 


f n = carrier frequency in the case of the carrier loop ("1.6xlO y Hz) 

u 7 

= P code rate in the case of the code loop ('1x10 Hz) 


In the s domain (Laplace transform), the pulse in jerk is given by 

-st, 

_ n _ , 

.4 


e(s) 


9 n-r->] 


The approach taken here will be to solve for the time response for 
a step in jerk (i.e., ignore the e" stl term) and then use superposition 
to get the pulse response by subtracting the step response at t= tj. The 
phase error, in the s domain, for the step in jerk is found to be 


A 2 - 
7 - 9 

4 >(s) ~ 9 


A 5 

3 A , In-* A — 

Z-! + ^ 0 * 7 es 

1 + As + ~ s 2 

F 


where A = 


r +1 

2W. 


The inverse transform is first found for the case of critical damping; 
i.e. , c = 1 


0 f 25 


’ 3 164 

B l 


B L t (l+exp(- | B L t)) - HI (l-exp(- | B L t))]u(t) 


The more interesting case is for 0.707, for which we find 

♦ (t) « & IV ’ J ( Uexp( " I B L t) cos jB L t)]u(t) 


^9 \ 9 [ B. t 


Since this is the response to a step in jerk, the pulse response is found 
from superposition by using 

^^pulse “ U (t-V 

The normalized second-order loop error response for both a step 
input of jerk and a Dulse input (10 msec wide) are shown plotted in Fig- 
ure 21. It can be seen from the figure that, for a step long enough in 
duration, the loop error will grow to the point where the loop breaks 

lock. For a pulse input, the loop error reaches a steady-state value 

3 

that is proportional to the inverse of B L . Thus, the desire for wider 
tracking bandwidths to accommodate dynamic tracking can be seen. 

The GPS carrier loop and code loop errors for the Shuttle pitch 
start jerk input are shown plotted in Figure 22. The code loop bandwidth 
of = 1 Hz has purposely been chosen to be somewhat narrower than the 

actual case to accomplish a worst-case evaluation. Even so, the code 

_2 

loop error is seen to be only 10 P code chips, or about 1 ft. This 
error is insignificant. The carrier loop bandwidth of 10 Hz is also 
somewhat narrow for purposes of worst-case evaluation. Even so, the 

_3 

carrier loop error is seen to be only about 10 radians. This error is 
totally insignificant. Thus, it is concluded that typical on-orbit 
maneuvers do not affect the GPS receiver tracking performance. 




Figure 21. Error Response for 2nd Order Phase Loc 




Figure 22. GPS Receiver Tracking Loop Response to Shuttle Transient Maneuver 
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